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METHOD FOR CHARACTERIZING AND SIMULATING A CHEMICAL 
MECHANICAL POLISHING PROCESS 

Cross-Reference to Related Application : 

This application is a continuation of copending International 
Application No. PCT/DE01/04903 , filed December 27, 2001, which 
designated the United States and was not published in English. 

Background of the Invention : 
Field of the Invention : 
The present invention relates to a method for characterizing 
and simulating a chemical mechanical polishing process, in 
which a substrate that is to be polished, in particular a 
semiconductor wafer, is pressed onto a polishing cloth and is 
rotated relative to the latter for a defined polishing time. 

Chemical mechanical polishing (CMP) is a method for 
20 planarizing or polishing substrates that is in widespread use 
in particular in semiconductor manufacturing. By way of 
example, planarized surfaces have the advantage that a 
subsequent exposure step can be carried out with a higher 
resolution, since the required depth of focus can be lower on 
25 account of the reduced surface topography. 
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In this context, the problem arises that different structure 
densities and spacings in the layout of a semiconductor chip 
influence the planarizing properties of the CMP process. 
Inappropriately selected process parameters then lead to a 
5 considerable fluctuation in the layer thickness (global 

topography) over the chip surface after the CMP process. On 
the other hand, an unfavorably selected circuit layout also 
leads to insufficient planarization . In this context, the 
inadequate planarization, on account of the associated layer 

10 thickness variations over the chip surface or the image field 
surface of a subsequent exposure step, has an adverse affect 
on the subsequent processes and therefore also on the product 
properties. In particular the process window of a subsequent 
lithography step is reduced in size on account of the reduced 

15 depth of focus. 

Hitherto, the process parameters to be set for the CMP process 
have generally been adapted specifically for each new layer to 
be polished on the semiconductor wafer and for almost every 

2 0 new product. For each CMP process there are numerous process 
parameters, such as the rotational speeds of polishing plate 
and substrate holder, the compressive force, the polishing 
time, the condition of the polishing cloth or the choice of 
polishing abrasive. Furthermore, the deposition thickness of 

25 the layer which is to be planarized has to be matched to the 
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planarization properties of the CMP process used and the 
structure densities and sizes of the chip layout. 

The optimum parameters are typically determined in a series of 
5 test gradings by trial and error. These tests entail a not 
inconsiderable time and cost outlay and also require a 
sufficient number of wafers of a new product layout to be 
available . 

10 Furthermore, it is difficult to measure the resultant global 

topography on the test wafers, and consequently in practice it 
is often only the less relevant local planarization properties 
that are analyzed . 

15 Summary of the Invention : 

It is accordingly an object of the invention to provide a 
method for characterizing and simulating a chemical mechanical 
polishing process that overcomes the above-mentioned 
disadvantages of the prior art methods of this general type, 

2 0 in which the CMP process can be characterized in such a manner 
that for a given product layout the process result can be 
predicted without carrying out tests on real layout 
substrates . 

25 The method according to the invention for characterizing and 
simulating a CMP process, in which a substrate which is to be 
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polished, in particular a semiconductor wafer, is pressed onto 
a polishing cloth and is rotated relative to the latter for a 
defined polishing time, includes the steps of: defining a set 
of process parameters, in particular compressive force and 
5 relative rotational speed between the substrate and the 

polishing cloth; preparing and characterizing a test substrate 
having test patterns with different structure densities at the 
defined process parameters; determining a set of model 
parameters for simulating the CMP process from the results of 

10 the characterization of the test substrate; determining layout 
parameters of the substrate which is to be polished; defining 
a profile of demands on the CMP process result for the 
substrate which is to be polished; and simulating the CMP 
process in order to determine the polishing time required to 

15 satisfy the profile of demands. 

The method according to the invention has the advantage that 
an experimental characterization only has to take place once 
for a specific set of process parameters, specifically on a 
20 test substrate that has test patterns with different structure 
densities. The results of the characterization of the test 
substrate are used to determine a set of model parameters with 
which the CMP process can then be simulated for any desired 
layout . 

25 
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Then, for a given layout layout, parameters which form input 
variables for the simulation are determined. The demands 
imposed on the process result, for example a certain 
approximation to the optimum achievable global step height, 
5 are also defined. By simulating the CMP process, it is then 
possible to determine the polishing time required for this 
layout from the generally applicable model parameters and the 
specific layout parameters without experimental test grading 
using the layout itself being required. 

10 

Therefore, it is possible to determine on a theoretical basis, 
without using product wafers, whether a selected layout can be 
polished in the desired way using a specific process. It is 
also possible to reach conclusions as to the CMP process 
15 window. Therefore, the result is a considerable saving on 
time and costs in the technological development of new 
products . 

The test patterns of the test substrate contain regions with 
20 high (up) areas and low (down) areas of a defined step height, 
for example isolated blocks or line patterns. The ratio of up 
areas to down areas determines the structure density, the 
limits of which are formed by a density of 0% (only down 
areas) and 100% (only up areas) . A preferred test substrate 
25 includes line patterns with a period (the width of the up and 
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down areas together) of 250 |im for structure densities of 4% 
to 72%. 



In one configuration of the method, the test substrate is 
5 characterized in an experimental polishing time grading in 

which the layer thickness development of the test patterns is 
measured in dependence on the polishing time. 

Preferably, the set of model parameters determined contains 
10 the abrasion rate, the hardness of the polishing cloth, and a 
characteristic filter length for determining effective 
structure densities. In this case, an effective structure 
density is obtained from the specific structure density of a 
layout by determining or forming a suitable mean over an area 
15 of a certain size. 

It is preferable for the mean to be formed by convolution of 
the specific structure density with a weighting function. The 
weighting function selected is expediently a two-dimensional 

20 Gaussian distribution, and the characteristic filter length is 
in this case the half -width value of the Gaussian curve. 
However, there are also other suitable weighting functions, 
for example square, cylindrical and elliptical weighting 
functions. According to current knowledge, the elliptical and 

25 Gaussian weighting functions have the minimum errors and are 
therefore preferably used. 
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The abrasion rate and the hardness are advantageously 
determined from the layer thickness development of a test 
pattern with a mean structure density. In this context, it is 
5 expedient for the abrasion rate to be determined from the 
pitch of the layer thickness development for long polishing 
times, and for the hardness of the polishing cloth to be 
determined from the speed at which the up and down areas of 
the test patterns reach the abrasion rate. The values for the 
10 abrasion rate and the hardness can, for example, be obtained 
by matching a local polishing model to the experimental 
results of a polishing time grading. 

The filter length is advantageously determined from the global 
15 step height of all the test patterns on the test substrate. 
In this case, the global step height is the difference in 
layer thickness between the maximum layer thickness measured 
value for all the up areas and the minimum layer thickness 
measured value for all the down areas. Since the global step 
2 0 height therefore represents a correlation over the surface of 
the entire layout, it is quite plausible that a significant 
global step height may remain even though the local steps have 
already been leveled by the polishing operation. However, it 
is the global step height over the image field area of a 
25 subsequent exposure step (for example 21 x 21 mm 2 ) that is 
crucial to the depth of focus of the exposure step. 
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In one configuration of the method, the layout parameters of 
the substrate used are the minimum and maximum effective 
structure density, p m i n and p maX / and the starting step height. 
5 The effective structure densities in turn result from the 

specific structure density of the layout by forming a suitable 
mean over an area of a certain size, characterized by the 
filter length. 

10 In a further configuration, a surface coverage with structures 
is determined for at least one region on the substrate, in 
order subsequently to use a cross-sectional profile of the 
corresponding structures to calculate a local structure 
density from the surface coverage and the cross-sectional 

15 profile of the structures. This is because the starting 

topography that is to be planar i zed by a CMP process is not 
determined by the layout directly, but rather is also 
determined by the preceding processes, such as for example an 
etching or deposition process. 

20 

In this context, account is taken of the fact that, by way of 
example, structures which have been etched in a preceding 
process or covered with a layer no longer have a box-shaped or 
rectangular profile, but rather on the one hand have an edge 
25 which is set back or projects with respect to its base and on 
the other hand also have angled or curved edges. Recesses or 
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angled edges for a given surface coverage lead to a reduction 
in the structure density compared to box-shaped structures of 
precise surface area and therefore also to a reduction in the 
amount of material to be removed, while projecting edges lead 
5 to an increase. The effective structure density is then 
calculated by forming the mean over the filter length. 

The simulation method therefore also takes account of the 
preceding process. For a given structure having a width and a 

10 height, it is possible to cite a cross-sectional profile for a 
specific known preceding process. To do this, it is possible 
to store corresponding measured data in tables in order for 
them then to be assigned to the structures of the existing 
surface coverage during the simulation, or alternatively it is 

15 possible to cite simplified geometric formulae which are 

applied to the corresponding profile of the structure below. 

To calculate the local structure density, a first volume is 
calculated by integration of the cross-sectional profile over 

20 the basic area of a structure and then the first volume is 
divided by a second volume, which is calculated from the 
product of the basic area of the structure and the starting 
height. Given a mathematically predetermined function of the 
cross-sectional profile, the integration can be carried out 

25 directly, or alternatively numerical integration is carried 
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out by use of nested intervals. The two integrals converge as 
the number of interval steps moves towards infinity. 



The profile of demands that has been defined is preferably 
5 given by a global step height to be achieved on the substrate 
after the CMP process has been carried out, since the global 
step height has a crucial influence on the depth of focus of a 
subsequent exposure step . 

10 In one configuration of the simulation method, the deposition 
thickness required to carry out the CMP process is determined 
in addition to the required polishing time in the simulation. 

The simulation preferably also determines the minimum global 
15 step height that can be achieved. This determination is based 
on the discovery that for sufficiently long polishing times 
the local steps have disappeared and the global step height 
only changes to a negligible extent. For the limit scenario 
of an infinitely long polishing time, the result is a residual 
20 global step height which is dependent only on the starting 

step height and on the minimum and maximum effective structure 
density which can be achieved in the layout which is to be 
polished. 



25 If the minimum achievable step height is determined, it is 

recommended for the global step height that is to be achieved 
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to be selected as a function of the minimum achievable global 
step height. By way of example, working on the basis of the 
starting step height, it is required to achieve 80%, 90% or 
95% of the difference between the starting step height and the 
5 minimum achievable global step height. A procedure of this 
type represents a compromise between being sufficiently close 
to optimum planarization and the demand for short polishing 
times . 

10 The invention also includes a method for the chemical 

mechanical polishing of a substrate, in particular of a 
semiconductor wafer, in which a CMP process is simulated as 
described, a layer which is to be planarized is deposited on a 
substrate and the substrate is polished for a polishing time 

15 determined from the simulation. As has been stated, it is not 
necessary to carry out a new experimental test grading for 
each new substrate layout. Rather, the results of an 
experimental characterization of a test substrate can be used 
for a wide range of product layouts. 

20 

In the polishing method, the CMP process is preferably 
simulated using a method that also provides the required 
deposition thickness as a simulation result. The layer that 
is to be planarized is then deposited in the required 
25 thickness before the polishing step. 
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Other features which are considered as characteristic for the 
invention are set forth in the appended claims. 

Although the invention is illustrated and described herein as 
5 embodied in a method for characterizing and simulating a 

chemical mechanical polishing process, it is nevertheless not 
intended to be limited to the details shown, since various 
modifications and structural changes may be made therein 
without departing from the spirit of the invention and within 
10 the scope and range of equivalents of the claims. 

The construction and method of operation of the invention, 
however, together with additional objects and advantages 
thereof will be best understood from the following description 
15 of specific embodiments when read in connection with the 
accompanying drawings . 

Brief Description of the Drawings : 

Fig. 1 is a flow diagram of a CMP simulation method according 
2 0 to the invention; 

Fig. 2 is a flow diagram illustrating a subroutine of the flow 
diagram shown in Fig. 1 in more detail; 
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Fig. 3 is a graph plotting a measured layer thickness in the 
up area and down area of a structure of average density and 
also a global step height as a function of polishing time; 

5 Fig. 4 is a graph plotting the measured global step height and 
the global step height obtained from the CMP simulation model 
as a function of the polishing time; 

Fig. 5 is a diagrammatic illustration relating to the 
10 definition of sizes used in a CMP polishing process; 

Fig. 6 is a cross-sectional profile of a substrate with 
structures on which an HDP deposition process has been carried 
out ; 

15 

Fig. 7A is a graph illustrating the layer thickness applied in 
an HDP process as a function of the lateral extent of a 
structure; 

2 0 Fig. 7B is a graph illustrating the layer thickness applied in 
an HDP process as an integration of the profile by nested 
intervals; 

Fig. 8A is a plan view of the surface coverage of two windows 
25 with structures before an HDP process; 
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Fig. 8B is a plan view of the surface coverage of two windows 
with structures after the HDP process; 

Fig. 9A is a graph showing a diagram as in Fig. 7 but for a 
5 conformal deposition process; and 

Fig. 9B is a graph showing a diagram as in Fig. 7 but for an 
etching process. 

10 Description of the Preferred Embodiments : 

Referring now to the figures of the drawing in detail and 
first, particularly, to Fig. 5 thereof, there is shown 
diagrammatically , to define the sizes used, a wafer 12 which 
is to be polished and a polishing cloth 18. The wafer 12 has 

15 a structure containing high up areas 14 and low down areas 16 
with a step height h 0 . On account of the rotational movements, 
a local relative speed v results between the wafer 12 and the 
polishing cloth 18 at any location. A compressive force F and 
a surface area of the wafer 12 can be used to determine a 

20 local abrasion rate in a known way using the Preston's 
equation. 

Fig. 1 shows a flow diagram of an exemplary embodiment of a 
chemical mechanical polishing (CMP) simulation method 100. In 
25 a first step 102, a relative speed of the wafer 12 and the 
polishing cloth 18 and the compressive force, for example a 
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relative rotational speed or table speed (TS) = 35 rpm 
(revolutions per minute) and a compressive force of 6 psi, are 
defined as process parameters of the process which is to be 
characterized. 

5 

In step 104, the selected process is completely characterized 
as a one-off. To do this, as illustrated in the flow diagram 
presented in Fig. 2, first a suitable test substrate is 
selected (reference numeral 2 02) . In the exemplary 

10 embodiment, the test substrate has test patterns containing 
isolated blocks and line patterns with different structure 
densities of 4% to 72%. All the structures of the test 
patterns have relatively large dimensions (> 10 \im) in order 
to allow simple optical examination of the structures and to 

15 enable their development to be evaluated as a function of the 
polishing time. 

The test substrate is characterized in step 204, the result 
obtained being the layer thickness development for various 
20 structure densities as well as the global step height as a 
function of the polishing time (reference numeral 206) . 

Then, in steps 206 to 214, the experimental values are 
reproduced by use of a local CMP model with a global density 
25 by matching model parameters abrasion rate K, polishing-cloth 
hardness E and filter length cO . 
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The abrasion rate K and the hardness of the polishing cloth E 
are determined from the layer thickness development of a test 
pattern of average structure density, as illustrated in Fig. 
5 3. 

Fig. 3 plots the measured layer thickness in the up area 
(reference numeral 3 02) and down area (reference numeral 3 04) 
of a structure of average density. It can be seen that 
10 substantially only the high, up area is abraded, while the 
abrasion rate in the down area is low. 

At slightly longer times, the down area is also abraded, and 
for relatively long polishing times the abrasion rates for the 
15 up and down areas converge (reference numeral 310) . The pitch 
of the layer thickness curves in the area 310 is then a 
measure of the abrasion rate K. 

The hardness E of the polishing cloth determines how quickly 
2 0 the up and down areas reach the abrasion rate. The precise 

values for K and E are determined by matching a local model to 
the results of the polishing time grading. The details of a 
local model of this type are described, for example, in the 
article titled "A CMP Model Combining Density And Time 
25 Dependencies'' by Taber H. Smith et al . , Proc . CMP-MIC, Santa- 
Clara, CA, Feb. 1999. 
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The filter length cO is obtained from the development of the 
global step height over the course of time. The global step 
height is in this case the layer thickness difference between 
5 the maximum layer thickness measured value of all the up areas 
and the minimum layer thickness measured value of all the down 
areas at each time, 

Stgiobai(t) = Max Up - Min Down . (1) 

As can be seen from the plot of the measured global step 
height 306 illustrated in Fig. 3, the global step height is 
still significant when the local step height, i.e. the 
difference between the layer thickness in the up area 
(reference numeral 3 02) and the layer thickness in the down 
area (reference numeral 3 04) has already virtually disappeared 
for a test structure of defined structure density. 

The CMP model is now matched to the profile of the global step 
20 height by obtaining an effective structure density p(x,y), 

which is likewise included in the model calculation, from the 
specific structure density p 0 (x,y) of the test substrate by 
convolution with a weighting function. 
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Each weighting function in this case has a characteristic 
filter length cO, which indicates the size of the area used to 
form the mean. In the exemplary embodiment, the weighting 
function selected is a two-dimensional Gaussian distribution 
5 with a half-width value cO . 

It has now been found that for given process parameters the 
global step height St g i oba i(t) which remains, given sufficiently 
long polishing times, is dependent only on the starting step 
10 height h 0 and on the minimum and maximum effective densities of 
the layout, in this case of the test substrate : 

St g iobal(t->«) = h 0 (p ma x-pmin) (2) 

15 Since p max and p m i n are dependent on cO, the filter length can 
be determined by comparing equation (2) with equation (1) for 
sufficiently long times . 

In the model calculation, therefore, the value of the filter 
20 length cO is a fit parameter which is iteratively adapted 
until the simulated data sufficiently match the data 
determined experimentally in the polishing time grading (steps 
208, 210, 212, 214) . 
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Fig. 4 shows the result of a CMP simulation after adjustment 
of the filter length cO . Fig. 4 illustrates the measured 
global step height 4 02 and the global step height 4 04 obtained 
from the model as a function of the polishing time. 

5 

At the end of the process characterization 104, the model 
parameters K, E and cO have been matched to the selected 
process conditions. The result is then a simulation model 
that can be applied to any desired product layout without 
10 further free parameters . 

Returning now to Fig. 1, in step 106 layout parameters are 
determined for specific application to a product layout. For 
this purpose, the minimum and maximum effective densities of 
15 the product layout and the starting step height are determined 
from the specific structure density of the product layout, 
which is known from measurements or from the design data, by 
use of the weighting function with the filter length cO . 

20 A simulation of the CMP process for the product layout using 

the previously determined values for K, E and cO then directly 
results in the local and global step heights as a function of 
the polishing time. 

25 As can be seen from the global step height plotted in Fig. 4, 
the global step height does not drop to zero over the course 
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of time, but rather, after a sufficiently long polishing time, 
tends toward its limit value given by equation (2) . There is 
therefore no point in continuing polishing for a very long 
time, since this lengthens the process time without 
5 significantly improving the process result. 



Therefore, in step 106 of the simulation method, a profile of 
demands imposed on the CMP process result is defined; 
satisfying the profile of demands results in that the 
10 polishing process can be ended. For this purpose, in the 

exemplary embodiment a variable o is determined, for example 
at a value of 0.95, indicating what proportion of the maximum 
achievable polishing result is sufficient for the specific 
polishing process. 

15 

This cessation condition then enables the CMP simulation to 
determine the polishing time t p i an required. This results from 
the equation 



20 Stgiobal (t plan ) " St glo bal ( t ->•) = (l-O) (h 0 -St g iobal ( t->- ) ) , 

i.e. for a = 0.95, the global step height is reduced by 95% of 
the maximum possible reduction from h 0 within the polishing 
time tpian* 
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Furthermore, a layer thickness S d0 wn which has been abraded in 
the down area with the lowest effective structure density at 
the time t p i an can be used to determine the deposition thickness 
A required to achieve this degree of planarization : 

5 

A = Sdown ( tpian, pmin) + ho 

Therefore, the material thickness which is to be applied, the 
required planarization time and the resulting global step 
10 height can be determined by the simulation without it being 
necessary to use real product wafers. 

In an alternative exemplary embodiment, to determine the 
effective structure density by subtraction or addition of 
15 critical structure sizes from the surface coverage p' (x,y) 

according to the chip layout which are characteristic of the 
preceding processes and for subsequent surface coverage 
determination, the density of the surface topography of the 
structures following the preceding process is determined. 

20 

In this case, the specific structure density during the 
deposition is defined as the ratio of volume to the product of 
a window area 400 of individual structures or of a field of 
structures under consideration and the maximum step height h 0 . 
25 In the case of precisely one structure, this corresponds to 
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the basic area of the structure. Since the filter length cO 
of a CMP process is approximately 1 mm, it is possible for the 
window areas 400 within which this surface coverage is 
determined to be selected to be small compared to the filter 
5 length cO but large compared to an individual structure. 

An exemplary embodiment considered here is an algorithm for 
determining the HDP deposition topography on a metal level. 
Fig. 6 shows a typical determined cross-sectional profile of a 
layer 302 deposited in this manner. The HDP deposition is 
used to fill trenches with a high aspect ratio. Structures 
with a lateral size below a defined dimension (on the left- 
hand size in Fig. 6) are grown over completely, with the 
result that flattened down areas 14' of a new surface 
topography are formed. More oxide is deposited on structures 
that are larger (on the right-hand side in Fig. 6) , so that up 
areas 14' that have been changed from the structure layout are 
formed, and a flank 15' is formed at their edges. The flank 
15' is characteristic of the HDP process. It changes with the 
process parameters of the HDP process. 

If the deposition height is plotted against the lateral 
structure size (Fig. 7A) , the result, in addition to the angle 
301 of the flank 15' , is two further characteristic lateral 
2 5 variables L min and L max - Lmin is half the lateral dimension below 
that a uniform deposition thickness grows over all the 
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structures of the structured metal layer. The thickness is 
the deposition height on an unstructured surface, reduced by 
the trench depth. Structures with a lateral extent of twice 
L ma x in turn have a constant deposition thickness grown over 
them and form a trapezoid (on the right-hand side of Fig. 6) . 
In this case, the height of the trapezoid is the deposition 
thickness on an unstructured surface. The structures between 
twice L min and L max are characterized in the profile by their 
pointed triangular shape (middle of Fig. 6) . The relationship 
between structure size and deposition thickness can in turn be 
defined by simulation or by SEN images and can be stored. 

When using numerical methods, the window area 4 00 is shifted 
over the layout and the surface coverage p' (x,y) therein is 
determined. As down areas 16', the surfaces associated with 
the regions L min do not make any contributions to the effective 
structure density, even though they contribute to the surface 
coverage. The areas of the edges 15' which are assigned to 
the regions L max are divided, by nested intervals, into a 
number n of intervals 305 each of known basic areas and are 
each provided with a mean value for the local structure height 
(Fig. 7B) . An inner region once again has a plateau, i.e. the 
up area 14' of height h 0 with respect to the down area. The 
product of the individual partial areas and the associated 
local structure heights results in the volume taken up by the 
material of the layer 302. This is set in a relationship with 
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respect to a volume that results from the product of the 
height h 0 times the window area 400. 



The result of the example HDP process is illustrated in Figs. 
5 8A and 8B. Fig. 8A shows two different surface coverages with 
the same structure densities in in each case one window area. 
The structures of the layout, i.e. the up areas 14, are 
illustrated in solid black. Fig. 8B accordingly in each case 
shows the remaining structure contributions after the 

10 topography of the HDP deposition process has been taken into 
account. The edges provided with an angle of inclination 301 
are placed beneath in graded grey shades in Fig. 8B in order 
to provide a plan view of the nested intervals. The result 
for the HDP process is not only a reduction in the structure 

15 densities compared to the layout, but also that this reduction 
is determined as a function of the structure extent or size, 
as can be seen from a comparison between the two windows 4 00 
in Fig. 8B. The finer structures in the layout (smaller up 
areas 14 and down areas 16) even provide exclusively down 

20 areas 16' following the HDP process. 

Figs. 9A and 9B show further examples of processes with a 
deposition height or the structure height plotted against the 
lateral structure extent, specifically for a conformal 
25 deposition process on structures with a low aspect ratio (Fig. 
9A) and an etching process (Fig. 9B) . The, for example, 
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experimentally determined variables L m i n and L max and also H 0 may 
in this case also adopt negative values i.e. by way of example 
may have the effect of increasing the topography compared to 
the structure from the layout. 

5 

The determination of the layout parameters p m i n and p max as the 
minimum and maximum values for the effective structure density 
is carried out after the mean has been formed for the specific 
structure density having the filter length c 0 as calculated 
10 from the cross-sectional profile and the surface coverage. 

Of course, it is also within the scope of the invention to 
select a different set of process parameters, to carry out the 
CMP simulation using this set of parameters and to compare the 
15 results with those obtained above in order to optimally adapt 
the process parameters to a given product layout. 
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